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ABSTRACT 

We report on the results of the first simultaneous X-ray (RXTE) and optical (SAAO) 
observations of the luminous low mass X-ray binary (LMXB) GX9+9 in 1999 August. 
The high-speed optical photometry revealed an orbital period of 4.1958 hr and con- 
firmed previous observations, but with greater precision. No X-ray modulation was 
found at the orbital period. On shorter timescales, a possible 1.4-hr variability was 
found in the optical light curves which might be related to the mHz quasi-periodic oscil- 
lations seen in other LMXBs. We do not find any significant X-ray/optical correlation 
in the light curves. In X-rays, the colour-colour diagram and hardness-intensity dia- 
gram indicate that the source shows characteristics of an atoll source in the upper ba- 
nana state, with a correlation between intensity and spectral hardness. Time-resolved 
X-ray spectroscopy suggests that two-component spectral models give a reasonable fit 
to the X-ray emission. Such models consist of a blackbody component which can be 
interpreted as the emission from an optically thick accretion disc or an optically thick 
boundary layer, and a hard Comptonized component for an extended corona. 

Key words: accretion, accretion discs - binaries: close - stars: individual (GX9+9) 
- X-rays: stars 



1 INTRODUCTION 

GX9+9 (4U 1728-16) is a persistent bright galactic bulge 
X-ray source, with a 4.19 ± 0.02 hr orbital period that was 
reported by Hertz and Wood (1988) from HEAO-1 data. 
The variations were roughly sinusoidal with an amplitude of 
3.8%. However, there was clearly substantial scatter (4.6% 
residual error) about this modulation. The source showed 
a rich variety of X-ray variability on timescales from min- 
utes to months and correlated X-ray spectra and fluxes (Ma- 
son et al. 1976; Parsignault & Grindlay 1978). Observations 
with EXOS AT show a spectrum and spectral variations very 
similar to GX3+1, GX13+1 and GX9+1 (Schulz, Hasinger 
& Triimper 1989) and so GX9+9 was classified as a high- 
luminosity atoll source (Hasinger & van der Klis 1989) . The 
X-ray spectra obtained by EXOSAT and Einstein are very 
similar and can be represented by a blackbody of 1.5 keV and 
a thermal bremsstrahlung component of ~ 8 keV (Christian 
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& Swank 1997). An ASCA observation revealed a similar 
kind of spectrum with a blackbody of ~ 1.3 keV plus a cut- 
off power law with photon index of ~ 1.2 and cutoff energy of 
~ 6 keV (Church & Balucihska- Church 2001). Line emission 
was also found from previous observations {Einstein: 0.77, 
0.82 and 1.06 keV; Christian & Swank 1997; ROSAT: 0.65 
and 0.97 keV; Schulz 1999) but a recent ASCA observation 
found none (Church & Balucinska-Church 2001). An RXTE 
observation of GX9+9 was carried out in 1996 but no kHz 
QPOs were found (Wijnands, van der Klis & van Paradijs 
1998). There is no report of X-ray bursts from this source. 

The optical counterpart (V ~ 17) is a 'Sco X-l-like' blue 
star with He II A4686 and C III/N III A4640 emission lines 
(Davidsen, Malina & Bowyer 1976; Charles, Thorstensen & 
Bowycr 1977). The only temporal study in the optical per- 
formed to date is that of Schaefer (1990) in which a ~ 20% 
modulation was present on the X-ray periodicity. There is no 
radio detection for GX 9+9 according to the recent Australia 
Telescope Compact Array observations at 4.8 and 8.64 GHz 
(Berendsen et al. 2000). 
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There have been no simultaneous optical/X-ray obser- 
vations of this source so far. At V = 17 it is ideal for studying 
outer disc variability simultaneously with RXTE observa- 
tions of the X-ray variability. In this paper, we report the 
results of the first simultaneous X-ray/optical observations 
of GX9+9. An observation log for both X-ray and optical 
observations is presented in while the temporal analysis 
of the data is in In ^\ time-resolved X-ray spectroscopy 
is presented and an archival study of X-ray modulation is 
reported in §5. The implications of our new data will be 
discussed in iJS| 



2 OBSERVATIONS 

2.1 SAAO Observations 

All optical observations were carried out with the UCT-CCD 
fast photometer (O'Donoghue 1995) at the Cassegrain focus 
of the 1.9-m telescope at the South African Astronomical 
Observatory (SAAO), Sutherland from 1999 August 11 to 
14. The observing conditions were good with typical seeing 
of ~ 1-2 arcsec and there was ~ 6 hr observation time each 
night. In general, the seeing during the first two nights was 
quite stable while the image quality of the remaining two 
nights was relatively poor due to the variable seeing and 
high wind conditions. In order to achieve high time reso- 
lution, no filter was employed (i.e. white light) and a 5 x 5 
on-chip prebinning was applied to maximise the time resolu- 
tion. The timing resolution of all the data is 2 s. In Tabled 
an observing log is presented together with the X-ray ob- 
servations. Debiasing and flat-fielding were performed with 
standard IRAF routines. Since the star fields are crowded, 
differential photometry was obtained via PSF fitting with 
IRAF. The light curves of all four nights are shown in Fig.0 
to 0] and have typical error bars on individual 2 s observa- 
tions of less than 2%. 

Apart from the high-speed photometry, we also ob- 
tained a set of B and V images on 1999 August 15. The expo- 
sure times on GX 9+9 were 60 s in the B-band and 40 s in the 
U-band. We calibrated the magnitudes using a photometric 
standard, giving B = 16.96 ± 0.05 and V = 16.79 ± 0.05. 

2.2 RXTE Observations 

GX 9+9 was observed with the Proportional Counter Array 
(PCA) instrument on board the RXTE (Bradt, Rothschild 
& Swank 1993) between 1999 August 11 and August 14 (see 
Table . The PCA consists of five identical Proportional 
Counter Units (PCUs) sensitive to X-rays with an energy 
range of 2-60 keV and a total effective area of ~ 6500 cm 2 . 
The PCA spectral resolution at 6 keV is approximately 18 
per cent, and the maximum timing resolution available is 
1 [is. Three PCUs were on during our observations. For most 
of the time, RXTE observed the source simultaneously with 
the optical telescope at SAAO. The typical on-source time 
for each night is about 10-14 ks (see Table0 . Apart from the 
standard configuration of the PCUs (Standardl and Stan- 
dard2 modes which are always available), we also employed 
event modes with maximum timing resolution of 125 fjs and 
64 spectral energy channels in order to maximise the timing 
and spectral resolution. All X-ray spectra presented here 



Table 1. A Journal of the X-ray/optical observations of GX9+9. 



Observatory Detector Date Time (UT) 









Start 


End 


SAAO 


1.9m + 


1999 August 11 


17:26 


23:08 




UCT CCD 


1999 August 12 


17:15 


23:08 






1999 August 13 


17:16 


23:00 






1999 August 14 


17:45 


23:05 


RXTE 


PCA 


1999 August 11 


17:20 


00:06 






1999 August 12 


17:28 


00:03 






1999 August 13 


17:27 


22:41 






1999 August 14 


17:26 


22:40 



have been corrected for background and dead-time with 
FTOOLS v4.2 and Epoch 4 PCA background model. 



3 TIMING ANALYSIS 
3.1 Optical 

It is very clear from Figs. to [I] that substantial variability 
is present in the optical light curves of GX9+9. In order 
to search for any modulation in the light curves, we used 
the Lomb-Scargle periodogram (Lomb 1976; Scargle 1982; 
hereafter LSP), a modification of the discrete Fourier trans- 
form which is generalised to the case of uneven spacing. We 
have also employed Phase Dispersion Minimisation (PDM; 
Stellingwerf 1978) to check the results from LSP, as PDM 
is more sensitive to non-sinusoidal modulation. Before per- 
forming period-searching, the light curves were normalised 
to fractional count rate in order to eliminate nightly varia- 
tions. All the light curves were then combined into a single 
dataset which was used to calculate the LSP. Fig. |S] shows 
the LSP of the combined light curve and there is a very 
strong peak at 4.1958 ± 0.0005 hr. Also present in the LSP 
are the aliases from this modulation. This is demonstrated 
by applying the LSP to a simulated light curve of a pure 
sinusoidal modulation at 4.1958 hr and with the same time 
sampling as the real data (see Fig. [5J. This shows how the 
LSP ought to appear if there were only a 4.1958 hr mod- 
ulation in the data. The principal features of the LSP are 
restored when comparing with the real data. A folded light 
curve on the 4.1958 hr modulation is shown in Fig. in 
which the phase zero is set at the time of the first data 
point (HJD 2451402.23164) and it is highly sinusoidal with 
semi-amplitude of 9.8%. Figure 6 shows the nightly folded 
light curves on the 4.1958 hr modulation and the sinusodial 
modulation is very stable. It is worth noting that the phase 
zero roughly corresponds to minimum light as determined 
from the data. A PDM analysis was also performed which 
confirmed these results. The LSP of several check stars with 
similar brightness as the target star were also calculated but 
none of these show any modulation at 4.1958 hr, with an 
upper limit of semi-amplitude ~ 0.17% . Such small scale 
variations are mainly due to the local atmospheric trans- 
parency variations commonly found at SAAO (e.g. Homer, 
Charles & O'Donoghue 1998). The combined results indi- 
cate that the highly sinusoidal and significant modulation 
at 4.1958 hr of GX9+9 represents the intrinsic behaviour of 
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Figure 1. Simultaneous X-ray and optical light curves of GX9+9 on 11/8/99. The X-ray light curve was divided into two energy bands: 
2-3.5 keV (upper panel) and 9.7-16 keV (middle panel); typical error bars for the X-ray data are also shown. The timing resolution is 
16 s for X-ray data and 2 s for optical data. For clarity, the optical light curve (lower panel) has been rebinned to 10 s. Note that a 4.2-hr 
orbital modulation is obvious in the optical light curve. Also shown underneath the source is a light curve with the 4.2-hr sinusoidal 
modulation subtracted. Optical light curves of a comparison star are also plotted. 



the source and is entirely consistent with the previously re- 
ported orbital period (Schaefer 1990) but to a greater degree 
of accuracy. 

For modulation on shorter timescales (~ 1000 s), we 
expect any sub-orbital modulation which is either periodic, 
or quasi-periodic will be evident in the LSP or PDM of the 
data, and should repeat from night to night. Therefore the 
LSP and PDM for each night were calculated in order to 
search for such variability. Before calculating both the LSP 
and PDM, the 4.19-hr orbital modulation was subtracted 
(assuming a sinusoid at this period; see Fig.^to^J. In Fig.|7| 
the nightly LSP is shown together with the 99% confidence 
level in both white noise (Gaussian noise) and red noise. It is 
apparent that in Fig. [7|those peaks between 0.5-1 hr _1 are 
below the red-noise level. From the log-log binned power 
density spectrum (PDS; insets of Fig. 0, both the actual 
and simulated PDS are similar, indicating the variability 
of the source on short timescales can be described as the 



combination of white and red noise. All the PDS level off 
at high frequency (~ 30 hr _1 ), suggesting that white noise 
becomes dominant between 30 and 900 hr _1 (the Nyquist 
frequency). There is also an indication that an excess ap- 
pears at ~ 1.4 hr from the log- log PDS. In determining 
the white-noise confidence level, we generate Gaussian noise 
datasets with the same time intervals and variance as the 
true data and then perform the LSP analysis on the result- 
ing datasets. The peak power in each periodogram (which 
must be purely due to noise) was then recorded. This was 
repeated 10,000 times for good statistics. However, the noise 
is not necessarily Gaussian, and can also be frequency de- 
pendent with a higher power at lower frequencies. By using 
the above method, strong peaks at the low frequencies might 
give misleading results. In order to take into account such a 
noise contribution in the data, we simulate the noise in the 
form of a power law, which is also known as red noise (e.g. 
Done et al. 1992) . Quantitatively, the power spectrum of the 
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Figure 2. Same as Fig. [I] but for 12/8/99. 



"red-noise" light curve will be given by (l/u>) a , where u is 
frequency and a is the slope index. Following an implemen- 
tation of the method of Timmer & Konig (1995), we generate 
the simulated light curve with a red-noise power spectrum, 
together with the same time sampling and variance as the 
original data. Finally, we combine the white-noise simulated 
light curve as mentioned above, with the red-noise simu- 
lated light curve, by scaling the data points in each such 
that the relative contribution of both white and red-noise 
components is as specified. Following the method used pre- 
viously to determine the white-noise confidence level, we de- 
rive the red-noise confidence level with the combined simu- 
lated light curve. Since red noise is frequency dependent, we 
need to compute the confidence level for a set of frequency 
bins. Hence it looks like a histogram rather than a horizontal 
line as for the white-noise level (see Fig. 0. Note that the 
peak at around 1.4 hr (0.7 hr _1 ) is just below the red-noise 
level in the first night yet it seems to repeat itself in other 
nights, demonstrating its quasi-periodic nature. Hence, the 
negative result from the LSP suggests that such variabili- 
ties are not strictly periodic, but are due to some form of 



quasi-periodic behaviour. The 1.4-hr variability may be an 
intrinsic property of the source and it may be associated 
with the mHz QPOs which were also found in other LMXBs 
(e.g. Chakrabarty et al. 2001). 



3.2 RXTE 

A similar temporal analysis was performed for the RXTE 
data. In order to study the variability in different energy 
bands, light curves from two energy ranges (2-3.5 keV and 
9.7-16 keV) were extracted for each night. From Fig. 
it is clear that flaring is stronger at higher energies. In soft 
X-rays, the count rate roughly remains at a constant level, 
while the flares in hard X-rays can be by a factor of ~ 2. 
The LSP and PDM were calculated from the combined 2-16 
keV band pass, as in the case of the optical data. The PDM 
method is particularly important for the X-ray data since 
the source changed its spectral state from time to time (see 
for details about the spectral behaviour of the source) 
and therefore any modulation would be expected to be non- 
sinusoidal. Fig.|H|shows the resulting LSP and PDM for the 
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Figure 3. Same as Fig. [I] but for 13/8/99. 



2-3.5 keV RXTE data. There is no distinct peak at the op- 
tical periodicity of 4.19 hr (marked) and hence the 4.19-hr 
X-ray modulation is below the detection threshold in our ob- 
servations. At the 3-cr confidence level, upper limits of 0.9% 
and 6.7% can be set for the semi-amplitude of any 4.19-hr 
variability in the 2-3.5 keV and 9.7-16 keV bands respec- 
tively. The 3.8% amplitude seen in HEAO-1 data is below 
the hard X-ray detection threshold and it is not implausi- 
ble that a weaker signal is present. Note that the sudden 
jump at about 0.6-0.9 hr -1 in the 99% red noise confidence 
level is due to the window function arising from the data 
gaps, typically seen in low-Earth orbit satellites. Since the 
red noise is frequency dependent, the confidence level will 
be different for different frequency bins and therefore such 
analysis also reflects the window function (if it exists) in the 
data. 

Variability on shorter timescales is also seen in the LSP 
but it is not significant. A period analysis for each night was 
also performed and the resulting LSPs are shown in Fig. |5] 
There is no significant peak at the 4.19-hr orbital period 
which confirms the result from the combined light curves. 



On shorter timescales, those marginally significant peaks do 
not repeat more than one night and they are also comparable 
to the observing windows of the spacecraft orbit, typically 
< 1 hr due to SAA passages. Hence they are not considered 
further. 

The X-ray PDS of GX9+9 was also calculated for each 
night in order to search for kHz QPOs. The PDS was de- 
rived from the 125 /is time resolution data over 100-4000 Hz. 
However, the analysis resulted in a negative finding and an 
upper limit (95% confidence) was set of 1.6% to the rms 
amplitude of QPOs with a FWHM of 150 Hz (in 2-60 keV). 



3.3 Comparison of Light Curves 

The simultaneous X-ray/optical light curves (Fig. ex ~ 
hibit complex variability and it is interesting to look for 
any correlation. We first applied cross-correlation analysis 
to compare the sinusoidal subtracted optical light curves 
with the X-ray light curves in different energy bands. We 
also divided the light curves into smaller intervals to study 
the cross- correlation. We repeated this analysis for all four 
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Figure 4. Same as Fig. [I] but for 14/8/99. 



observations. There is no significant correlation indicating 
that the correlated response is very weak. Fig. I1UI shows a 
close up of the light curves from Fig. in which there is a 
prominent X-ray flare in the second RXTE visit. However, 
there is no immediate response in the optical. Similar uncor- 
rected variability can also be found from other light curves. 
We have tried to use the transfer function technique as em- 
ployed by Kong et al. (2000) to investigate the correlations, 
but the optical light curves could not be reproduced as a 
convolution of the X-ray light curve with a Gaussian trans- 
fer function. Hence this suggests that the optical emission 
does not significantly respond to X-ray variability. We note, 
however, that the presence of uncorrelated events does not 
imply the absence of correlated ones. 



4 X-RAY SPECTROSCOPY 

Colour-colour (CD) and hardness-intensity diagrams (HID) 
were derived from the same four energy bands as used by 
Wijnands, van der Klis & van Paradijs (1998) to facilitate 
comparison: 2-3.5, 3.5-6.4, 6.4-9.7 and 9.7-16 keV. The soft 



and hard colours are defined as the count rate ratio between 
3.5-6.4 keV and 2-3.5 keV and between 9.7-16 keV and 6.4- 
9.7 keV, respectively. The total intensity for the HID was 
taken as the sum of the total count rates (1 PCU) in the 
same energy bands as were used in calculating the colours. 
All the data were binned to a time resolution of 64 sec- 
onds. The CD and HID of all the X-ray data are shown in 
Fig. EH Both the CD and HID of GX9+9 suggest that the 
source was on the banana branch (Hasinger & van der Klis 
1989) during the RXTE observations. Both diagrams show 
very clear correlation and this is typical for GX 9+9 (see e.g. 
Schulz, Hasinger & Triimper 1989; Wijnands, van der Klis 
& van der Paradijs 1998). 

As there is clear indication of correlation between in- 
tensity and colour ratio from the CD and HID, correlated 
spectral parameters and intensity are expected, therefore 24 
approximately half-hour summed spectra were derived from 
the light curves (see Table |5J for further spectral fitting. 
Based on the spectral analysis of archival Crab Nebula and 
pulsar data (e.g. Dove et al. 1998; Wilms et al. 1999; Bloser 
et al. 2000a; Barret et al. 2000; Kong 2000), we restrict our 
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Figure 6. Nightly folded optical light curve of GX9+9 on a period of 4.1958 hr. Two cycles are shown for clarity. Phase zero is arbitrarily 
set at the time of the first data point (HJD 2451402.23164). The relatively smaller amplitude of the August 13, 1999 data is due to the 
bad seeing near the end of the run (see also Figure 3). 
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Figure 7. LSP of GX9+9 after subtracting the 4.19-hr sinusoidal modulation. The dotted horizontal line is the 99% confidence level 
assuming white noise, while the dotted histogram is the 99% confidence level due to red noise. The 4.19-hr orbital period is marked by 
a vertical dashed line. Insets: the log-log binned power density spectra (with error bars), with the average of 200 simulated white- noise 
plus red- noise models over-plotted (thick line). Note that the power spectra level off at high frequency due to the large contribution of 
white noise. 



spectral analysis to the energy range of 2.5-20 keV, outside 
which the calibration is not reliable. In addition, a 1% sys- 
tematic error was added to account for the uncertainty of 
the calibration before the spectral fitting. 

Previous attempts to fit the X-ray energy spectrum 
of GX9+9 show that it cannot be fitted well with a 
single-component model (either power law, blackbody, cut- 
off power law or thermal bremsstrahlung; e.g. Schulz 
1999). These single-component models were fitted first to 
our RXTE data but the results were similarly negative 
with unacceptable (> 2). Two-component models were 
then applied to the spectra. Combinations of blackbody 
plus cutoff power law (BB+CPL), blackbody plus ther- 
mal bremsstrahlung (BB+TB), and blackbody plus Comp- 
tonization (BB+CompST; Sunyaev & Titarchuk 1980 and 
BB+CompTT; Titarchuk 1994) result in a significant im- 
provement to the spectral fits (see Tables [3])^] and Fig. I12H . 
All models included a column density Nh, which was kept 



fixed at 2.1 x 10 21 cm -2 (e.g. Church & Balucihska-Church 
2001) for all the fits since it tended to vary over a wide range 
if left as a free parameter. We also note that from an optical 
spectroscopic observation (Shahbaz et al. 1996), the inter- 
stellar absorption features suggest a reddening of E(B — 
V) = 0.1 — 0.4. Hence, the column density can be derived by 
using the conversion of iV w /10 22 cm- 2 = 0. 179,4 v (Predehl 
& Schmitt 1995) and A v = 3.1E(B - V) (Schild 1977), sug- 
gesting Nh is in the range of 0.6 — 2.2 x 10 21 cm -2 . The total 
flux of GX9+9 varies between 6.8-8.9 xlO" 9 erg cm" 2 s~\ 
corresponding to 2 — 2.7 x 10 37 D 2 ergs" 1 where Ds is in 
units of 5 kpc (e.g. Church & Balucihska-Church 2001) over 
the four days of our observations. The best-fit parameters 
of all the models for the 24 spectra are given in Tables HJ- 
|S] The effective radius of the blackbody emitting surface 
Rbb and the ratio of the blackbody flux to the total flux 
in the 2-20 keV band, Fbb/Ftotai are also included. For the 
BB+CompST and BB+CompTT models, the Comptonizing 
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Figure 8. LSP and PDM analysis of the combined RXTE data in the 2—3.5 keV band. The horizontal dotted line is the 99% confidence 
level from white noise (WN), while the dotted histogram is the 99% confidence level due to red noise (RN). Note that the sudden jump 
in the red noise confidence level is due to the window function of the data. 
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Figure 9. LSP of the nightly RXTE data in 2—3.5 keV. The dotted horizontal line is the 99% confidence level from white noise, while 
the vertical dashed line is the 4.19-hr orbital period. 



y-parameter, defined as y = 4fcT e r 2 / 'm e c 2 was derived. In 
Fig. 1121 the spectral parameters of all the models are plot- 
ted as a function of the total (2-20 keV) flux. For all the 
models, the blackbody temperature is about 1.4-2 keV and 
is always correlated with the total flux. As the total flux in- 
creases, the blackbody component increases from ~ 10% to 
40% while the other component remains roughly constant. 
Hence, the increase of the blackbody emission may also be 
responsible to the flaring activities. Such positive correlation 
also reflects the correlation between the intensity level and 



the hardness ratios described in the HID. The blackbody 
radius Rbb drops slightly from ~ 2.5 km to ~ 2 km when 
the total flux increases. 

For the BB+CPL model, the photon index is roughly 
at ~ 1.6 with no indication of any correlation with the total 
flux. Meanwhile, the cutoff energy E cll t varies from 5.5 keV 
to 7 keV. Such results indicate that the contribution of the 
CPL component remains roughly constant as the total flux 
increases, while the BB component becomes the major driv- 
ing force for the spectral evolution of the source. For the 
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Figure 10. A close up of the simultaneous X-ray and optical (sinusoidal subtracted) light curves of GX9+9 on 11/8/99. The timing 
resolution is 1.875 s for X-ray data and 2 s for optical data. The count rate for both datasets is normalised to fractional count rate. 



BB+TB model, the thermal bremsstrahlung temperature 
varies from 5 keV to 6 keV and is correlated with the total 
flux. 

For the BB+CompST, the optical depth r is moderate 
at about 9-12 and might be roughly anti-correlated with the 
blackbody contribution. The electron temperature is very 
steady at ~ 3 keV although small variations can be seen 
and are also anti-correlated with the optical depth. The 
BB+CompTT results are very similar to BB+CompST but 
with larger error bars in r and kT e . The additional parame- 
ter in CompTT, the seed photon temperature (To), is steady 
at ~ 0.5 keV. Also r in CompTT is slightly smaller than that 
of the CompST while kT e is larger than that of the CompST. 

We also investigate the effect of replacing the BB com- 
ponent with a multi-colour disc blackbody model (DBB; 
Mitsuda et al. 1984) in all models. All the fits were accept- 
able but resulted in a large colour temperature for the inner 



disc (kTin > 3 keV) and an unfeasibly small inner disc ra- 
dius (Rin = 1 km). The error bars for the other component 
(i.e. CPL, CompST and CompTT) were also large compared 
to the BB models. Therefore the DBB model seems to be 
unrealistic. 



5 ARCHIVAL X-RAY OBSERVATIONS 

If the 4.19-hr modulation is orbital in nature, it should be 
very stable (as in the optical). In order to study the X-ray 
modulation on this timescale, light curves from the public 
data archives were obtained and only observations longer 
than 15 ks ( > 4.2 hr) were considered. TableHJshows all the 
relevant observations with the date and duration of the ob- 
servations. An ~ 17.4-hr nearly uninterrupted EXOSAT ob- 
servation and an ~ 13-hr A SCA observation were extracted 
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Figure 11. The colour-colour diagram (a) and hardness-intensity diagram (b) of all four nights data of GX9+9. Each data point is an 
accumulation of 64 s. 



Table 2. Spectral fits of GX9+9 with the BB+CPL model. 



Number 


Start (UT) 


End (UT) 


pa 


(keV) 




^bb 
(keV) 


r tot 


F bb /Ff ot 




R bb 
(km) 


xl 

(34 dof) 


1 


11 


17:22:23 


17:55:59 


1.43 ±0.07 


5.85 ±0.39 


2 


04 ± 0.04 


8.3 





23 


1 


61 ±0.06 


0.68 


2 


11 


18:29:03 


18:54:19 


1.59 ± 0.10 


6.41 ±0.58 


1 


69 ± 0.02 


7.61 





21 


2 


12 ± 0.05 


0.88 


3 


11 


19:02:23 


19:35:59 


1.44 ±0.08 


5.60 ±0.38 


1 


88 ± 0.04 


7.82 





20 


1 


69 ± 0.07 


0.67 


4 


11 


20:25:43 


21:15:59 


1.65 ± 0.08 


6.13 ±0.41 


1 


59 ± 0.02 


7.08 





18 


2 


14 ±0.05 


1.08 


5 


11 


21:15:59 


22:39:19 


1.66 ± 0.07 


6.48 ± 0.39 


1 


58 ± 0.01 


7.11 





21 


2 


31 ± 0.03 


0.60 


6 


11 


23:12:23 


23:45:59 


1.71 ± 0.10 


6.68 ±0.59 


1 


63 ± 0.02 


7.41 





21 


2 


27 ±0.06 


0.91 


7 


12 


14:28:15 


18:01:51 


1.65 ± 0.10 


6.10 ±0.53 


1 


57 ±0.02 


6.89 





16 


2 


07 ±0.05 


0.77 


8 


12 


18:34:55 


19:25:11 


1.64 ±0.07 


6.06 ±0.35 


1 


58 ± 0.02 


6.95 





16 


2 


00 ±0.05 


0.68 


9 


12 


19:58:15 


20:31:51 


1.61 ± 0.10 


5.93 ±0.49 


1 


58 ± 0.03 


6.90 





14 


1 


90 ± 0.07 


0.88 


10 


12 


20:31:52 


21:05:11 


1.60 ±0.09 


5.97 ±0.43 


1 


59 ± 0.02 


7.12 





18 


2 


13 ± 0.05 


0.91 


11 


12 


21:38:15 


22:45:11 


1.53 ± 0.07 


5.56 ±0.29 


1 


55 ± 0.02 


6.76 





13 


1 


86 ± 0.05 


0.79 


12 


13 


17:26:23 


17:59:59 


1.64 ±0.10 


6.09 ±0.52 


1 


56 ± 0.02 


7.19 





19 


2 


34 ± 0.06 


1.01 


13 


13 


18:16:23 


18:49:59 


1.42 ± 0.09 


5.48 ±0.57 


2 


26 ± 0.07 


8.77 





28 


1 


47 ±0.09 


0.87 


14 


13 


18:50:00 


19:39:59 


1.53 ± 0.07 


6.31 ±0.41 


2 


05 ± 0.04 


8.42 





25 


1 


65 ± 0.06 


0.5 


15 


13 


19:56:23 


20:29:59 


1.59 ± 0.10 


6.11 ±0.51 


1 


68 ± 0.03 


7.44 





20 


2 


07 ±0.07 


0.86 


16 


13 


20:29:59 


20:54:59 


1.61 ± 0.10 


6.46 ±0.55 


1 


73 ± 0.02 


7.80 





24 


2 


19 ± 0.05 


1.21 


17 


13 


21:36:23 


22:09:59 


1.52 ± 0.08 


6.16 ±0.45 


1 


97 ± 0.04 


8.16 





24 


1 


74 ± 0.07 


0.86 


18 


13 


22:09:59 


23:00:00 


1.77 ±0.10 


7.56 ± 0.75 


1 


72 ± 0.02 


7.56 





24 


2 


21 ± 0.05 


1.07 


19 


14 


17:25:35 


17:59:11 


1.60 ±0.09 


7.07 ±0.63 


2 


05 ± 0.04 


8.93 





29 


1 


85 ± 0.07 


0.82 


20 


14 


18:23:55 


18:50:51 


1.55 ± 0.08 


6.76 ±0.60 


2 


14 ± 0.04 


8.75 





27 


1 


60 ±0.06 


0.78 


21 


14 


18:57:15 


19:39:11 


1.60 ±0.09 


6.39 ±0.51 


1 


72 ± 0.02 


7.67 





22 


2 


10 ±0.05 


0.64 


22 


14 


19:55:35 


21:02:31 


1.73 ± 0.07 


6.76 ± 0.42 


1 


60 ±0.01 


7.15 





21 


2 


31 ± 0.03 


0.7 


23 


14 


21:35:35 


22:09:11 


1.61 ± 0.10 


6.55 ±0.56 


1 


58 ± 0.02 


7.29 





19 


2 


28 ± 0.06 


0.9 


24 


14 


22:17:15 


22:42:31 


1.58 ± 0.10 


6.82 ±0.74 


1 


89 ± 0.05 


8.00 





21 


1 


75 ± 0.09 


0.93 



a Photon index 
b Cutoff energy 
c Blackbody temperature 

d Total flux in 2-20 keV (10~ 9 erg cm" 2 s" 1 ) 
e Ratio of blackbody flux to total flux, 2-20 keV 
f Blackbody radius, assuming a distance of 5 kpc 
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Figure 12. Variation of spectral parameters for the BB+PL, BB+TB, BB+CompST and BB+CompTT models with the total flux in 
2-20 keV. 
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Table 3. Spectral fits of GX9+9 with the BB+TB model. 



Number 




kT bb 


kT t% 


Ftot 


Fbb/Ftot 




Rbb 








(keV) 


(keV) 








(km) 


(35 dof) 


\ 


o 


07 4- n 05 


c rc; -|_ n 1 A 


8.27 


21 


4 


AO, 4. o 07 
^±y zc u.u/ 


0.69 


2 


4 


f.a 4- n 05 


AQ 4- OQ 
u.^ty m u.uy 


7.55 


15 


4 


77 i 0.10 


0.92 


3 


4 


no _|_ o 05 


r: _|_ n i i 
U . OU ZC U . ± -L 


7.79 


0.17 


4 


4. o 08 
u^ zc u.uo 


0.66 


4 


4 


54 -1- 05 


t; HQ + Of! 
u.uy m u.uu 


7.00 


0.11 


4 


75 zb 0.10 


1.25 


5 


4 


52 + 04 

Uii _1_ U.Urt 


u.ijy _i_ u.uu 


7.10 


0.13 


4 


as 4- o 10 


0.72 


5 


4 


□! U.UJ 


U . ^U U . U t 


7.33 


13 


4 


76 ± 0.10 


1.06 


7 


4 


p;q _i_ n ofi 


U . U f u.u/ 


6.82 


10 


4 


71 ± 13 


0.90 


$ 


4 


_|_ n OR 


r: o 1 ? 4- D^i 
u.uo m u.uu 


6.88 


0.10 


4 


64 ± 0.13 


0.80 


9 


4 


uo m u.uo 


r: r\e: _)_ o 07 
u . uu m u.u/ 


6.84 


0.09 


4 


51 ± 0.15 


0.92 


10 


4 


uu m u.uu 


c i o _i_ o Ofi 

U . J. ^ IC U . UU 


7.06 


11 


4 


77 ± 0.11 


1.06 


44 


4 


51 i 0.07 


c: 01 4- Ci^ 
U . Ul u . uu 


6.71 


08 


4 


54 ± 0.14 


0.84 


12 


1 


53 ± 0.05 


5.06 ± 0.07 


7.12 


0.12 


1 


93 ±0.13 


1.11 


13 


2 


31 ± 0.04 


5.20 ±0.26 


8.73 


0.28 


1 


41 ±0.05 


0.83 


14 


2 


12 ± 0.05 


5.62 ±0.15 


8.36 


0.21 


1 


42 ±0.07 


0.58 


15 


1 


68 ± 0.05 


5.28 ±0.08 


7.38 


0.14 


1 


73 ±0.10 


0.95 


16 


1 


74 ± 0.04 


4.48 ± 0.09 


7.73 


0.18 


1 


87 ±0.08 


1.30 


17 


2 


02 ± 0.05 


5.58 ±0.13 


8.11 


0.20 


1 


51 ±0.07 


0.91 


18 


1 


72 ± 0.05 


5.58 ±0.09 


7.47 


0.15 


1 


71 ±0.10 


1.42 


19 


2 


13 ± 0.04 


5.87 ±0.18 


8.87 


0.25 


1 


57 ±0.06 


0.92 


20 


2 


24 ± 0.05 


5.79 ±0.21 


8.70 


0.24 


1 


37 ±0.06 


0.86 


21 


1 


72 ± 0.05 


5.45 ± 0.08 


7.61 


0.16 


1 


77 ±0.10 


0.75 


22 


1 


56 ± 0.04 


5.25 ±0.06 


7.07 


0.12 


1 


83 ± 0.09 


0.92 


2.3 


1 


53 ± 0.05 


5.51 ±0.07 


7.23 


0.13 


1 


99 ±0.13 


1.03 


24 


1 


99 ± 0.07 


5.75 ±0.15 


7.94 


0.16 


1 


36 ±0.10 


0.95 



a Thermal bremsstrahlung temperature 



from the archival database. In addition, we also analysed our 
BeppoSAX data taken in early 2000. Both the LSP and PDM 
were applied in the same way to search for any periodicity. 
Fig. ESI shows the results from the LSP for the ASCA GIS, 
EXOSAT ME and BeppoSAX LECS data. Clearly no pe- 
riod at 4.19 hr is detected. In addition, variability on shorter 
timescales is not seen in the LSP, which confirms our RXTE 
results described above. 



6 DISCUSSION 

6.1 Orbital Modulation 

We have observed the luminous LMXB GX 9±9 simultane- 
ously with X-ray and optical telescopes. In the optical, the 
source was seen at B — 16.96 ± 0.05 and V = 16.79 ± 0.05, 
which is very consistent with previous observations (Schae- 
fer 1990). The high quality optical light curves of GX9+9 
presented here confirm the 4.19-hr periodicity found previ- 
ously (Schaefer 1990). The stability of this period also shows 
that it is the orbital period of the system. Unlike in the op- 
tical, no significant X-ray modulation on the 4.19 hr period 
was detected in the four-night RXTE observations and other 
archival X-ray observations (see §5). This is somewhat in- 
consistent with the previous HEAO-1 results (Hertz & Wood 
1988), in which a 3.8% modulation was found. In fact, the 
HEAO-1 data showed large scatter (4.6% residual error) 
and the light curve was based on 74 data points, typically 
~ 1.9 hr apart and spread over 5.91 days. Their detected 
4.19 hr X-ray modulation is only marginally significant at 
the 99% confidence level. Comparing the two observations 




0.4 0.6 0.8 1 

Frequency (1/hr) 



Figure 13. LSP of the ASCA, EXOSAT and BeppoSAX obser- 
vations of GX 9+9 (see Tablel6l. The dotted horizontal line is the 
99% confidence level, while the vertical dashed line is the 4.19-hr 
orbital period. No significant peak at the orbital period is seen in 
either plot. Note that the peaks in the BeppoSAX LSP are due 
to the window function arising from the data gaps, typically seen 
in low-Earth orbit satellites. 
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Table 4. Spectral fits of GX 9+9 with the BB+CompST model. 



Number 




kT bb 
(keV) 


XT? 
(keV) 




y c 


Ftot 


Fbbl Ftot 




Rbb 

(km) 


(34 dof) 


1 


I 


79 4- n ri3 
. < ^ zn u.uo 


9 nc; 4_ n 1 n 

i.yj Z\Z U. J-U 


11 99 4- n 59 


2.90 


8.35 


0.17 


I 


.97 ± 0.07 


0.69 


2 


I 


c;n 4- n ri4 


3 nn 4- n 1 3 

O.UU HI U. J-O 


1 n 36 -1- 53 

J-U.OU ZH KJ.OO 


2.50 


7.64 


0.21 


2 


fiq + n 14 
. uy zn u. j-^i 


0.90 


3 


I 


t;K 4 (1 (It 
- ju zn u.uo 


9 S3 4- n OQ 
i.oo in u.uy 


11 99 4- n 49 


2.78 


7.86 


0.17 


2 


31 xnnq 

. O -L IC U. U3 


0.69 


4 


I 




3.11 ± 0.11 


Q 40 -1- 34 


2.15 


7.13 


0.22 


2 


27 4 (1 (18 


1.17 


5 


I 


4K _|_ n 09 

■ rtU _1_ U . 


3 1 fi 4- n 09 

J.1U _1_ \J.\JiJ 


Q 4S + 28 

C/.t:0 _1_ U.iiO 


2.22 


7.16 


0.24 


2 


-yo _i_ u.uo 


0.62 


(] 


I 


cn -1- fl 09 


q 99 _i_ n 1 q 

. u. to 


q 97 4. tO 


2.16 


7.46 


0.23 


2 


.78 i 0.07 


0.96 




I 


-^i^ m u.uo 


q n« _i_ n 14 

o.uu u.-L'i 


Q _|_ 4S 

y . uo m u . 10 


2.18 


6.93 


0.20 


2 


.79 i 0.12 


0.88 




I 


.^±0 □! U . Ui 


q 1 9 4- n OQ 


q Q7 4. 9Q 
j.o/ nr. u . ^y 


2.14 


7.00 


0.20 


2 


ai 4. ns 
.01 m u.uo 


0.75 


9 


I 


4n -i- n ri9 


■3 09-1-0 10 
o.u^ nr. u.tu 


Q JK _|_ ^4 

y . 1 nr. u . o'i 


2.24 


6.94 


0.19 


2 


. Ol ' nz u. uo 


0.94 


10 


I 


40 J- n 09 


q on 4- in 

o.uu m u.tu 


07 4. t7 
y . / u . / 


2.28 


7.16 


0.23 


3 


or -I- n HQ 

.uu m u.uy 


1.13 


11 


I 


Q7 -1- n 09 

■ O I m u . u^ 


^ no 4- ns 

o.uu u.uo 


q 77 4. 9Q 
y . 1 / nr. u . ^y 


2.24 


6.81 


0.17 


2 


QO 4- n OR 

-OU □! U.UO 


0.91 


12 


1 


.43 ±0.02 


3.08 ±0.13 


9.51 ±0.42 


2.18 


7.23 


0.22 


2 


.99 ±0.09 


1.08 


13 


1 


.96 ±0.02 


2.81 ±0.11 


12.18 ±0.68 


3.60 


8.79 


0.10 


1 


.20 ±0.02 


0.87 


14 


1 


.74 ±0.04 


2.89 ±0.09 


11.54 ±0.48 


3.00 


8.44 


0.19 


2 


.04 ±0.10 


0.58 


15 


1 


.50 ±0.03 


3.03 ±0.13 


9.97 ±0.48 


2.35 


7.49 


0.22 


2 


.77 ±0.11 


0.96 


16 


1 


.56 ±0.03 


3.03 ±0.13 


10.23 ±0.52 


2.47 


7.83 


0.24 


2 


.70 ±0.10 


1.32 


17 


1 


.72 ±0.03 


3.08 ±0.14 


10.43 ±0.55 


2.62 


8.21 


0.20 


2 


.10 ±0.07 


0.91 


18 


1 


.58 ±0.03 


3.18 ±0.15 


9.73 ±0.48 


2.35 


7.59 


0.27 


2 


.75 ± 0.10 


1.35 


19 


1 


.86 ±0.02 


3.20 ±0.13 


10.41 ±0.48 


2.71 


8.97 


0.23 


2 


.01 ± 0.04 


0.86 


20 


1 


.87 ±0.03 


3.06 ±0.11 


11.15 ±0.49 


2.97 


8.78 


0.22 


1 


.93 ± 0.06 


0.83 


21 


1 


.52 ±0.04 


2.96 ±0.12 


10.46 ±0.50 


2.53 


7.71 


0.22 


2 


.74 ±0.15 


0.80 


22 


1 


.47 ±0.02 


3.12 ±0.09 


9.53 ±0.30 


2.22 


7.18 


0.23 


2 


.83 ±0.07 


0.82 


23 


1 


.47 ±0.02 


3.38 ±0.16 


9.12 ±0.43 


2.20 


7.34 


0.21 


2 


.76 ±0.07 


1.00 


21 


1 


.65 ±0.05 


3.16 ±0.20 


10.27 ±0.77 


2.61 


8.04 


0.18 


2 


.14 ±0.13 


0.93 



a Electron temperature of the Comptonizing cloud 

b Optical depth of the Comptonizing cloud 

c Comptonization parameter y = 4kT e r 2 /m e c 2 



from HEAO-1, the 4.19-hr X-ray modulation is only found 
in one of them. Since the source often has flaring and short- 
term variation (< 1 hr) especially in hard X-rays (> 10 keV; 
see Fig. such activities might cause a false alarm in 

determining the period and account for the large errors if 
the time resolution is poor. Hence, these results should be 
treated with caution. The RXTE data, and other archival ob- 
servations reported here provide much higher time resolution 
and energy-selected light curves to study the 4.19-hr modu- 
lation of this source in greater detail. Both LSP and PDM 
analysis show no indication of variability on this timescale. 
However, it is possible that an X-ray periodicity might be 
detectable from time to time as a result of changes in the ac- 
cretion geometry (such as the height of the disc rim) , and it 
should be noted that the X-ray period was discovered prior 
to the optical modulation. Consequently, the agreement be- 
tween the two periodicities means that the original X-ray 
modulation should not be discarded. 

The lack of X-ray modulation and eclipses or absorption 
dips on the orbital period is common amongst LMXBs and 
indicates that the inclination angle of the system must be 
low (< 70°). Hertz and Wood (1988) assume that a 1.4M Q 
neutron star is orbiting a low-mass main sequence secondary 
which is filling its Roche lobe and use the lack of eclipses to 
show that the inclination of GX 9±9 is less than ~ 63° . We 
also note that the observed semi-amplitude (9.8%) of the op- 
tical modulation suggests an orbital inclination of less than 
70° (van Paradijs & McClintock 1995). Although the RXTE 
data presented here cannot confirm the X-ray modulation 



found by Hertz and Wood (1988), the result is still consis- 
tent with the hypothesis that GX9±9 is a low inclination 
system. In order to account for the small amplitude X-ray 
variability on the orbital period, Hertz & Wood (1988) sug- 
gested that the disc does not hide the central neutron star 
(since it is at low inclination angle) and instead we are actu- 
ally observing directly the X-ray emission from the vicinity 
of the central source. Therefore Hertz & Wood (1988) pro- 
posed a variable covering of an accretion disc corona by the 
edge of the disc. Alternatively, Schaefer (1990) suggested 
that the X-ray modulation is caused by the changing visi- 
bility of X-ray flux reflected off the bright spot although it 
is more likely that this mechanism could explain the optical 
modulation. 



6.2 The Origin of the Optical Modulation 

The observed orbital period of GX9±9 provides a valu- 
able tool for understanding the system geometry. Assuming 
a canonical value for the mass of a neutron star (Mx = 
1.4 Mo), then from Kepler's third law, 

a = 1.471 i? (l + g) 1/3 (1) 

where a is the separation between the two stars, q is the 
mass ratio, i.e., MijMx and M2 is the mass of the compan- 
ion star. As the companion fills its Roche lobe, the radius of 
the companion is given by (Paczyhski 1971) 
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Table 5. Spectral fits of GX9+9 with the BB+CompTT model. 



Number 
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T 


y 
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F bb /F tot 




Rbb 

(km) 


(33 dof) 


\ 


I 


Qi _L n no 
- <y J. ZC U.UO 


n kk _|_ a no 
u.oo zc u.u^ 


A 1 O 4- f| QO 


7 nn _|_ 1 c,4 
t . yu zc j. . o'-t 


2.00 


8.19 


0.34 


2 


on 4. n f|7 

-ZrU ZC U.U J 


0.62 


2 


I 


c\A _|_ n nc 

• U^-t ZC U.LW 


fin 4(1 0^ 
u.ou zc u.uo 


Q CI _|_ fl CC 

o. o J. in u.uo 


Q QQ _|_ 1 Ofl 

O.OO ZC J- . ZrU 


2.04 


7.51 


0.30 


2 


.67 ± 0.16 


0.86 


3 


I 


7« _l n nc 

, f U ZC U.UJ 


n r:o _i_ n no 
u.uo zc u.u^ 


Q Cfi _|_ f| CI 
O. UU ZC U. O J. 


8.81 ± 1.37 


2.16 
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2 
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4 


I 
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2 
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5 
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6 
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^0 4- 09 
u.uu zc u.u^ 


t:. -Lt: ZC U. i U 
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2 
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7 


I 
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g 
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2 


.88 ± 0.18 


0.99 
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2 
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8.64 


0.34 


1 
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14 
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8.30 


0.31 


2 


.14 ±0.07 


0.50 
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0.31 


2 
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0.78 


16 


1 


.72 ±0.03 
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5.00 ±1.53 


6.50 ±1.69 


1.65 


7.68 


0.36 
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.71 ± 0.10 


1.14 


17 
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.90 ±0.03 
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5.76 ±2.42 
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3.06 


2 


.33 ± 0.07 


0.78 


18 
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.96 ±0.03 


0.53 ±0.02 
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0.75 
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7.56 


0.33 
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0.58 
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.59 ±0.03 


0.49 ± 0.02 
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7.53 ±0.76 


1.83 


7.04 


0.30 


2 


.77 ±0.10 


0.70 


23 


1 


.64 ± 0.04 


0.54 ± 0.02 


5.62 ± 1.90 


5.98 ± 1.66 


1.57 


7.17 


0.32 


2 


.73 ± 0.13 


0.79 


24 


1 


.75 ± 0.07 
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3.54 ±0.51 


9.45 ± 1.50 
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7.91 


0.25 


2 


.22 ± 0.17 
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a Temperature of seed photons 



Table 6. Archival X-ray observations of GX9±9. 



Date 


Satellite 


Energy 


Time span 


On-source time 


Remark 






Range (keV) 


ks / hr 


ks / hr 




1983 Sept 27 


EX OS AT 


0.05-2 


62.6 / 17.4 


51.7 / 14.4 


Nearly uninterrupted 


1994 March 25 


ASCA 


0.5-10 


47.2 / 13 


23.7 / 6.6 




2000 April 8 


BeppoSAX 


0.1-10 


100 / 27.8 


23.6 / 6.4 





shadow of the disc; and (iii) the bright spot in the edge of 
the accretion disc. However, we do not have unambiguous 
phasing information for this system and it is impossible to 
determine the variation of these components with respect to 
orbital phase. 

Haswell et al. (2000) proposed that GX9±9 is a pos- 
sible persistent superhumper based on the fact that o<0.3 
(based on its short period and compact object mass of at 
least 1.4Mq) and also the similarity of the fractional ampli- 
tude of the optical modulation to that of superhumps seen 
in SXTs (O'Donoghue & Charles 1996). Unlike in the case 
of cataclysmic variables, in which superhumps are caused by 
a tidally-driven modulation of the disc's viscous dissipation 
which varies on the beat between the orbital and disc preces- 
sion period, the optical emission of LMXBs is dominated by 
X-ray reprocessing. Hence a different mechanism is required 
to account for the superhump period. A possible solution is 
to associate the superhumps with varying X-ray reprocessing 
by changing disc area (Haswell et al. 2000). In other words, 
the disc surface area varies on the superhump period, re- 
sulting in a modulation in the optical emission. If this is the 



/ \ 1/3 

R 2 = 0.462a f j-L- j (2) 

From the main-sequence mass-radius relation (Patter- 
son 1984) 

(\ 0.88 
=1.3MR eq °- 8S (3) 

Solving the above three equations simultaneously, we 
obtain M 2 = 0.4 M , R Q = 0.45, q = 0.29 and a = 1.6 R e . 
This corresponds to an early M dwarf (M0 V-M5 V) with a 
surface temperature of ~ 3400 K and an absolute magnitude 
My ~ 10. From the measured V magnitude of 16.79, this 
gives Mv = 3 (assuming a distance of 5 kpc and Av = 0.3). 
Hence, the observed optical flux is substantially larger than 
that from the M star itself and it is clear that the optical flux 
must be due to X-ray reprocessing in the accretion disc and 
illuminated face of the companion star. There are at least 
three possible origins of the optical flux and its variation 
associated with the orbital period: (i) the accretion disc; (ii) 
the X-ray heating of the companion star which is not in the 
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case, we would expect a superhump period which is slightly 
(a few per cent) longer than the orbital period. However, 
with the only significant periodicity being that found in the 
optical (here and earlier work), the lack of an X-ray or spec- 
troscopic modulation still leaves open the possibility that 
GX9+9 is a persistent superhumper as the nature of this 
variation cannot be defined. However, if the X-ray period is 
true, we can rule out the superhump interpretation. 

6.3 Accretion Geometry 

The X-ray spectrum of GX9+9 has rarely been studied in 
detail even though it belongs to one of the few bright X-ray 
sources known before Uhuru (Bradt et al. 1968). In general, 
previous studies (e.g. Christian & Swank 1997; Church & 
Balucihska-Church 2001) and this work suggest that two- 
component models are more appropriate in fitting the energy 
spectrum. Such models consist of a blackbody component 
plus a Comptonized component (e.g. CPL, TB, CompST 
and CompTT). Similar behaviour is also found in other high- 
luminosity LMXBs (White, Stella & Parmar 1988). 

The physical interpretation of the origin of such multi- 
component models is still poorly understood although ob- 
servations have been made for many LMXBs (e.g., White, 
Stella & Parmar 1988; Barret et al. 2000; Bloser et al. 2000a; 
Bloser et al. 2000b). It is widely accepted that the hard com- 
ponent (i.e. the Comptonized component) is due to Comp- 
tonization of soft photons in a hot corona, while the source 
of the soft photons which lead to a Compton up-scattering 
and the soft thermal component are still the subject of de- 
bate (e.g. White, Stella & Parmar 1988). There are several 
possibilities to explain the origin of the soft thermal com- 
ponent. The most straightforward solution is that this soft 
component originates from an optically thick accretion disc 
or an optically thick boundary layer where the accreted mat- 
ter intersects the neutron star surface. However, Sunyaev & 
Shakura (1986) showed that the luminosity of the boundary 
layer should be at least equal to that of the disc if relativis- 
tic effects are included, whereas we have found that the BB 
contribution ranges only from ~ 10 — 40%. In fact, many 
authors have found similar results for other sources (White, 
Stella & Parmar 1988; Barret et al. 2000; Bloser et al. 2000a; 
Bloser et al. 2000b). Therefore Guainazzi et al. (1998) and 
Barret et al. (2000) suggest that the weak thermal compo- 
nent could be due to an optically thick accretion disc, while 
the other Comptonized component arises from an optically 
thin boundary layer. However, this picture is not suitable to 
explain the similarities in spectral behaviour between neu- 
tron stars and black holes since a black hole does not have a 
physical surface (see Barret et al. 2000). Moreover, the neg- 
ative result from the DBB model also suggests that the soft 
component might not be from the inner accretion disc, and 
such a model is not a good description of the RXTE data 
presented here. Therefore the Comptonized component is 
unlikely to be due to the boundary layer and such a disc 
corona model can likely be ruled out. 

Among several models mentioned here, the present data 
are more consistent with a picture in which the soft com- 
ponent is produced by an optically thick boundary layer 
while the hard component is from an extended hot corona. 
As discussed above, such a model cannot explain the small 
contribution of the blackbody emission. Hence some mech- 



anism should suppress the boundary layer emission relative 
to that of the corona and disc. One possibility is that the 
rotation period of the neutron star is close to equilibrium 
with the Keplerian period of the inner edge of the accretion 
disc (White, Stella & Parmar 1988) such that little energy is 
given up in the boundary layer. Another possibility is that 
part of the boundary layer emission is blocked by the in- 
ner disc corona or thickened disc (Lamb 1986; van der Klis 
et al. 1987). Such mechanisms have been discussed for the 
X-ray spectrum of 4U 1916-05 (Church et al. 1998; Bloser 
et al. 2000a). In particular, Bloser et al. (2000a) have shown 
that by fitting a BB+CPL model, the blackbody radius of 
4U 1916-05 decreases slightly as the accretion rate increases. 
In Fig. 1121 all four models show a decreasing trend for 
as the total flux increases. This may indicate that the in- 
ner disc can become thicker as the accretion increases and 
thereby obscure a larger region of the boundary layer. In 
addition, the derived Rbb here is consistent with ~ 2-3 km, 
implying that the blackbody component can come from only 
a limited region of the neutron star surface such as an op- 
tically thick boundary layer. It is also interesting to note 
that GX9+1 (White et al. 1986), 4U 1735-444 (Seon et al. 
1997) and 4U 1820-30 (Bloser et al. 2000b), which are all 
luminous atoll sources, show similar blackbody radii in the 
persistent spectrum. Spectral analysis of these sources also 
suggests that the blackbody component comes from an op- 
tically thick boundary layer. 

Applying the above accretion geometry to describe 
GX 9+9, one would expect to detect an iron K emission line 
at 6.4-6.7 keV which is a common feature of the spectra of 
many luminous LMXBs (e.g. White et al. 1986; Barret et al. 
2000), as a result of the accretion disc corona (see Fabian, 
Guilbert & Ross 1982). However, observations have shown 
that such emission is not present in GX 9+9 (see Church & 
Balucihska-Church 2001). Interestingly, GX9+1 is also well- 
known for its lack of iron K emission (White et al. 1986). A 
reflection model (pexrav in XSPEC) was also used to fit the 
spectra so as to search for a reflection component from the 
disc, but the derived reflection factor was consistent with 
zero and thus there is no evidence of reflection. One possi- 
ble reason for this is the relatively low corona temperature 
such that the hard photons from the corona will tend to be 
photoelectrically absorbed by the matter in the disc, rather 
than Compton scattered. The lack of the line emission and 
reflection component indicates that a significant fraction of 
the disc cannot be seen, either by reflection or irradiation. 

Regarding a physical interpretation of the spectral mod- 
els, a thermal bremsstrahlung interpretation for the hard 
component seems to be unrealistic since the derived ra- 
dius of the plasma cloud is too far from the neutron star 
(~ 10 9 cm) to be heated to several keV by the neutron 
star (White et al. 1986). On the other hand, although the 
CompTT model is an improved version of CompST which 
includes relativistic effects and extends to the cases of both 
optically thin and thick plasmas, the larger error bars of op- 
tical depth and electron temperature compared to CompST 
make the CompST model more suitable to describe the 
present RXTE data. Therefore, we will refer to the CompST 
model when discussing Comptonization below. For the CPL 
model, it approximates an unsaturated Comptonized spec- 
trum (y — 4kT e r 2 /m e c 2 ~ 1) of cool photons up-scattered 
in a hot electron plasma and its results may not be di- 
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rectly interpretable in terms of a physical description. There- 
fore the CompST model can give a more complete descrip- 
tion of the source spectrum. In fact, several luminous X-ray 
sources (such as 4U 1636-536, 4U 1820-30, ScoX-1, ScoX- 
2, CygX-2, GX9+1 and GX 17+2) have shown that the 
CompST and CPL models provide equally good fits to the 
source spectrum (Christian & Swank 1997). It is also in- 
teresting to note that the derived kT e and r for ScoX- 
1, 4U 1735-444, 4U 1916-05, 4U 1820-30, CygX-2, GX9+1, 
GX17+2 and KS 1731-260 are very similar to the results 
presented here (see Table [7J. 

The high-quality RXTE data presented here can be 
summarised with the following physical model: in the up- 
per banana state, the luminosity is relatively high and the 
blackbody radiation comes from an optically thick boundary 
layer which is partly obscured by the inner accretion disc, 
producing ~ 10-40% blackbody flux from a small equato- 
rial region of the neutron star. Increases in the blackbody 
flux are also responsible for the flaring activity. Similar be- 
haviour is also found in ScoX-1 (White, Peacock & Taylor 
1985) and ScoX-2 (White et al. 1986). A hot electron cloud 
(T ps 3.5 x 10 7 K) in the form of an accretion disc corona 
surrounds the neutron star and at such high luminosity, the 
blackbody photons from the boundary layer can heat up the 
corona to maintain its temperature (see Fig. I12H at about 
a constant level in spite of the Compton cooling. The rela- 
tively strong blackbody radiation produces more extensive 
ionisation and hence increases the optical depth r of the 
corona, reducing the radiation drag that generates QPOs 
in the framework of the sonic-point model (Miller, Lamb & 
Psaltis 1998), and thereby accounting for the lack of kHz 
QPO in GX9+9. It is also known that when other atoll 
sources (e.g. 4U 1636-53: Wijnands et al. 1997; 4U 1820- 
30: Smale, Zhang & White 1997) are in their upper banana 
branch, kHz QPOs do not occur. 



7 CONCLUSIONS 

Simultaneous RXTE and optical high-speed observations 
were performed to study the atoll source GX 9+9. The high- 
quality optical light curves from high-speed CCD photome- 
try during a four-night observation confirm, but with much 
greater precision, the 4.19-hr orbital period (Schaefer 1990). 
However, our RXTE observations cannot confirm the X-ray 
modulation at the orbital period (Hertz & Wood 1988). By 
using archival ASCA and EXOSAT data, and also our recent 
BeppoSAX observation, no X-ray modulation could be re- 
vealed. GX9+9 also shows variability on timescales ~ 1.4 hr 
in the optical light curves which could be related to the mHz 
QPO as seen in other LMXBs. We do not find any significant 
X-ray/optical correlation in the light curves. Observations 
indicate that the inclination angle of GX 9+9 cannot be too 
high (< 70°). The observed optical flux is mainly due to 
X-ray irradiation of the outer accretion disc and the repro- 
cessing of a bright spot in the disc and X-ray heating of the 
companion can contribute the observed modulation. The CD 
and HID show that GX 9+9 was in the upper banana state 
during observations. Time-resolved X-ray spectroscopy was 
also performed and several two-component spectral mod- 
els (BB+CPL, BB+TB, BB+CompST and CompTT) can 
provide a good fit to the energy spectra. A two-component 



spectral model is always required to fit the X-ray emission 
from several luminous LMXBs (see Table [7|. In general, 
such models consist of a blackbody component (soft com- 
ponent) and a Comptonized component (hard component). 
The blackbody component is often interpreted as the con- 
tribution from either an optically thick accretion disc or an 
optically thick boundary layer (interaction between the in- 
ner accretion disc and the neutron star's surface). On the 
other hand, the hard component may be from a hot extended 
corona, optically thin boundary layer. The likely interpreta- 
tion for GX 9+9 can be understood in terms of an optically 
thick boundary layer (soft component) which is partly ob- 
scured by the inner disc or thickened disc, while the hard 
Comptonized component is from an extended hot corona. 
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Table 7. BB+CompST spectral fits of other luminous LMXBs. 



Source kT^ kT e r Reference 

(kcV) (kcV) 



ScoX-1 


1.6 


3 


10.5 


White, Peacock & Taylor 1985 


4U 1735-444 


1.45-1.90 


3.05-3.26 


11.93-13.99 


Scon ct al. 1997 


4U 1916-05" 


1.49-2.03 


3.53-8.85 


5.78-13.46 


Bloscr ct al. 2000a 


4U 1820-30 a 


1.30-2.36 


2.75-7.94 


6.00-15.21 


Bloser et al. 2000b 


CygX-2 


1.13 


3.75 


8.4 


Christian & Swank 1997 


GX9+1 


1.49 


3.00 


10.9 


White, Stella & Parmar 1988 


GX17+2 


1.18 


3.00 


13.0 


White, Stella & Parmar 1988 


KS 1731-260 


1.00-1.30 


2.60-2.80 


10.40-12.60 


Barret ct al. 2000 


GX9+9 


1.37-1.96 


2.81-3.38 


9.12-12.18 


this work 



a Use CompTT instead of CompST for the Comptonized component 



Fabian A. C, Guilbert P. W., Ross R. R., 1982, MNRAS, 
199, 1045. 

Guainazzi M., Parmar A. N., Segreto A., Stella L., Dal Fi- 
ume D., Oosterbroek T., 1998, A&A, 339, 802. 

Haswell C. A., King A. R., Murray J. R., Charles P. A., 
2000, MNRAS, 321, 475. 

Hertz P., Wood K. S., 1988, ApJ, 331, 764. 

Homer L., Charles P. A., O'Donoghue D., 1998, MNRAS, 
298, 497. 

Kong A. K. H., Homer L., Kuulkers E., Charles P. A., Smalc 
A. P., 2000, MNRAS, 311, 405. 

Kong A. K. H., 2000, PhD thesis, University of Oxford. 

Lamb F. K., 1986. In: The Evolution of X- Ray Binaries, p. 
151, eds Triimper J., Lewin W. H. C, Brinkman W., 
Dordrecht: Reidel. 

Mason K. O., Charles P. A., White N. E., Culhane J. L., 
Sanford P. W., Strong K. T., 1976, MNRAS, 177, 
513. 

Miller M. C, Lamb F. K., Psaltis D., 1998, ApJ, 508, 791. 

Mitsuda K., Inoue H., Koyama K. et al., 1984, PASJ, 36, 
741. 

O'Donoghue D., 1995. Baltic Astronomy, 4, 519. 

O'Donoghue, D. & Charles, P. A. 1996, MNRAS, 282, 191. 

Paczyriski B., 1971, ARA&A, 9, 183. 

Parsignault D. R., Grindlay J. E., 1978, ApJ, 225, 970. 

Patterson J., 1984, ApJSS, 54, 443. 

Predehl P., Schmitt J. H. M. M., 1995, A&A, 293, 889. 

Schaefer B. E., 1990, ApJ, 354, 720. 

Schild R. E., 1977, AJ, 82, 337. 

Schulz N. S., Hasinger G., Triimper J., 1989, A&A, 225, 48. 



Schulz N. S., 1999, ApJ, 511, 304. 

Seon K.-L, Min K.-W., Yoshida K., Makino F., Lewin 
W. H. G., van der Klis M., van Paradijs J., 1997, 
ApJ, 479, 398. 

Shahbaz T., Smale A. P., Charles P. A., van Paradijs J., 
Hassall B. J. M., Callanan P., 1996, MNRAS, 282, 
1437. 

Smale A. P., Zhang W., White N. E., 1997, ApJL, 483, 119. 

Stellingwerf R. F., 1978, ApJ, 224, 953. 

Sunyaev R. A., Shakura N. I., 1986, Sov. Astron. Lett., 12, 
117. 

Sunyaev R. A., Titarchuk L. G., 1980, A&A, 86, 121. 
Timmer J., KSnig M., 1995, A&A, 300, 707. 
Titarchuk L., 1994. ApJ, 434, 570. 

van der Klis M. K., Stella L., White N., Jensen F., Parmar 
A. N., 1987, ApJ, 316, 411. 

van Paradijs J., McClintock J. E., 1995, In: X-ray Binaries, 
p. 58, eds Lewin W. H. G., van Paradijs J., van den 
Heuvel E. P. J., Cambridge Univ. Press, Cambridge. 

Vrtilek S. D., Penninx W., Raymond J. C, Verbunt F., Hertz 
P., Wood K., Lewin W. H. G., Mitsuda K., 1991, ApJ, 
376, 278. 

White N. E., Peacock A., Hasinger G., Mason K. O., Manzo 
G., Taylor B. G., Branduardi-Raymont G., 1986, 
MNRAS, 218, 129. 

White N. E., Peacock A., Taylor B. G., 1985, ApJ, 296, 475. 

White N. E., Stella L., Parmar A. N., 1988, ApJ, 324, 363. 

Wijnands R. A. D., van der Klis M., van Paradijs J., Lewin 
W. H. G., Lamb F. K., Vaughan B., Kuulkers E., 
1997, ApJL, 479, 141. 

Wijnands R., van der Klis M., van Paradijs J., 1998. In: 

IAU Sumposium 188 "The Hot Universe", p. 370, 



Simultaneous X-ray/optical observations of GX 9+9 



eds Koyama K., Kitamoto S., Itoh M., Dordrecht: 
Kluwer Academic. 

Wilms J., Nowak M. A., Dove J. B., Fender R. P., di Matteo 
T., 1999, ApJ, 522, 460. 



